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£h ' Abstract 

We first improve the previous results on the masses of the doublets (0 + , 1 + ) and 
(1 + , 2 + ) with QCD sum rules up to the order of 0(1/ mo) in HQET. With the mass 
matrix of the two 1 + states, we obtain the mixing angle for the two states and 
£0 • construct the interpolating currents which couple to only one of the two eigen states 

up to the order 0(1 /tuq). Using these results we derive the light-cone sum rules for 
| the amplitutes of the pionic decay of B\ and B^ up to the order 0(1/ mo) in HQET. 

£S| ■ The 0(l/m,Q) terms in these amplitutes are found to be small. When applied to D\ 

and L>2 the theoretical results are in good agreement with experimental data though 
\ the theoretical uncertainty is sizable for the D system. 



PACS number: 12.39.Hg, 13.25.Hw, 13.25.Ft, 12.38.Lg 



1 Introduction 

Important progress has been achieved in understanding the properties of heavy mesons 
composed of a heavy quark and a light quark with the development of the heavy quark 
effective theory (HQET) [|TJ. HQET provides a systematic expansion of the heavy hadron 
spectra and transition amplitude in terms of I/tuq, where tuq is the heavy quark mass. Of 
course one has to employ some specific nonperturbative methods to arrive at the detailed 
predictions. Among the various nonpeturbative methods, QCD sum rules is useful to 
extract the low-lying hadron properties pi. The propertes of the ground state heavy 
mesons have been studied with the QCD sum rules in HQET in 0, ||. In 0] the mass of 
the lowest excited heavy meson doublets (2 + , 1 + ) and (1 + , + ) were studied with QCD sum 
rules in the heavy quark effective theory (HQET) up to the order of 0(1/ mq). The widths 
for pionic decays of the lowest two excited doublets (2 + , 1 + ) and (1 + , + ) is calculated with 
conventional QCD sum rules in the leading order of HQET in ||. Properties of excited 
heavy mesons were investigated also with QCD sum rules in full QCD in 0. Deacy of 
(1 + , + ) were investigated with light-cone sum rules || in 0. The pionic couplings between 
the lowest three doublets (2 + ,l + ), (l + ,0 + ) and (1~,0~) of heavy mesons are investigated 
in HQET using the light-cone sum rules [|TUJ. In ]TT| the pionic couplings of the heavy 



baryons are calculated with the same technique. 
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One difficult problem encountered in studying the decay widths of excited heavy mesons 
with QCD sum rules is the following. Except for the lowest states 0~, 1~, the spectra 
contains a pair of states for any spin-parity J p . An interpolating current with definite 
spin-parity in general couples with both two states. This parity is especially serious for the 
two 1 + states. All theoretical calculations indicate that they are close in mass values but 
quite different in magnitude of their decay widths. One of the two 1 + states is a narrow 
resonance decaying mainly by emitting a D wave pion, while the other one is a very wide 
resonance decaying mainly by emitting a S wave pion. An interpolating current used for 
the narrow 1 + state with a small coupling to the other 1 + state may cause sizable error in 
the result of calculation. Only in the tjiq — > oo limit is there a conserved quantum number 
je, the angular momentum of the light component, which can be used to differentiate the 
two states. Therefore, HQET has the important and unique advantage for this purpose. 
Corrections for finite tuq can be calculated order by order in HQET. The interpolating 
currents which couple to only one of the two J p states have been given in and used in 

In this work we shall calculate the decay widths of the doublet (l + ,2 + ) to the first 
order in 1/mg using the light-cone sum rule approach. As the necessary steps we present 
the results from sum rules in HQET for masses of these excited states and their coupling 
constants to the interpolating currents at the leading and next-to leading order in section 
2 and section 3 respectivly. A large part of contents of these sections have been obtained 
in Q.The new things include the addition of the leading 0(a s ) perturbation term for the 
chromo-magnetic splitting £ and the results for corrections to coupling constants at the 
order ©(I/toq). We also improve the numerical analysis of the sum rules at the leading 
order. The continuum threshold u c and the working region for the Borel parameter T 
determined from the stability criterum of the sum rules in the leading order are used in 
the sum rules for £ term and the kinetic energy term K, in the corrections to the masses. 
Solving the mass matrix of the two distinct 1 + states, we arrive at the mixing angle for 
these two states. Using this angle the interpolating currents which couple to only one of 
the 1 + states up to the order 1/rriQ are constructed in section (|. In section | and | we 
employ the constructed currents and the light-cone QCD sum rules (LCQSR) || in HQET 
to calculate the pionic decay amplitudes of Bi and B\ up to the order 0{\/mq) taking 
into account the mixing of two 1 + states. The results are summarized and discussed in 
section |S]. 



2 Two-point QCD sum rules 

The proper interpolating current J^pj" 3 for the state with the quantum number j, P, ji 
in HQET was given in ||. These currents were proved to satisfy the following conditions 

{()./;:;;;;:'(()) ././"../;) = fp^p^^f 1 ^ , (i) 
<o|t (•/;:,;;:'(.'■)•/;/,:;''(())) io> = s^s^-iy s ■ ■ 

x f dt5(x-vt)U Pjh (x) (2) 
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in the mq — > oo limit, where 7] ai '" a ^ is the polarization tensor for the spin j state, v 
is the velocity of the heavy quark, = g a/3 — v a v 13 is the transverse metric tensor, S 
denotes symmetrizing the indices and subtracting the trace terms separately in the sets 
(«x • • ■ oij) and (Pi ■ ■ ■ Pj), fpj t and Hpj e are a constant and a function of x respectively 
which depend only on P and jg. Because of equations ([I]) and (H), the sum rule in HQET 
for decay widths derived from a correlator containing such currents receives no contribution 
from the unwanted states with the same spin-parity as the states under consideration in 
the rriQ — ► oo. Starting from the calculations in the leading order, the decay amplitudes 
for finite tuq can be calculated unambiguously order by order in the I/ttiq expansion in 
HQET. 

We shall confine our study to the doublets (0 + , 1 + ) and (l + ,2 + ) here. There are two 
possible choices for currents creating + and 1 + of the doublet (0 + , 1 + ), either 

4,+, 2 — h v q , (3) 

4r +>2 = ^ Kin* , (4) 

or 

j'0,+,2 = ~j= K(-i) p t q , (5) 

j'i + ,2 = 4| h v i 5 i?H) Ptq ■ (6) 

Similarly, there are two possible choices for the currents creating 1 + and 2 + of the doublet 
(1+2+). One is 

4,+,i = ^| hvl 5 (-i) (V? - \lt pt) q , (7) 

4,+vr = )fl~ h « ij Y- (V + irvr - \ gr 2 pt) ? . ® 

Another choice is obtained by adding a factor — i p t to (0) and @. Note that, without 
the last term in the bracket in (|7|) the current would couple also to the 1 + state in the 
doublet (0 + , 1 + ) even in the limit of infinite itiq. 

For the doublet (0 + ,l + ), when the currents J'o,+,i/2, ^'i,+,i/2 in ©, © are used the 
sum rule (same for the two states) after the Borel transformation is found to be 

f%*e- 2R/T = ^ "V^du; -^rnl (qq) . (9) 

The corresponding formula when the current Jo,+,2 and Ji,+,2 m © and (^) are used instead 
of </'o,+,2 and J'i,+,2 is the following 

^, 1/2 e- 2A/T = ^ jf w V^du + \ (qq) - ^ mg (ft) . (10) 
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When the currents (|7]) and (@) are used the sum rule for the (1 + , 2 + ) doublet is found to 
be 

Here m\ (qq) = (qga^G^q). These sum rules have been obtained in H]. In the derivations, 
only terms of the lowest order in perturbation and operators of dimension less than six 
have been included. 

For the QCD parameters entering the theoretical expressions, we use the following 
standard values 

(qq) = -(0.24 GeV) 3 , 
(a s GG) = 0.038 GeV 4 , 

m 2 = 0.8 GeV 2 . (12) 

Here we redo the numerical analysis of the sum rules with the following criteria. The high- 
order power corrections is less than 30% of the pertuabation term without the cuttoff uj c 
and the contribution of the pole term is larger than 40% of the continuum contribution 
given by the perturbation integeral in the region u > u c . u c is constrained to be smaller 
than 3.2 Gev which is roughly the expected position of the pole of the radial excited state 
with the same quantum numbers. Subjected to these criteria we arrive at the stability 
region of the sum rules u c = 2.8 - 3.2 GeV, T = 0.43 - 0.73 GeV for (0), u c = 2.8 - 3.2 
GeV, T = 0.82 - 1.1 GeV for and u c = 2.8 - 3.2 GeV, T = 0.54 - 0.73 GeV for flu]). 
The results for A for three cases are the following 

A(l/2 + ) = 1.15 ±0.10 GeV, (13) 
for the doublet (0 + , 1 + ) when the currents (|3|), (f|) without the derivative is used. 

A(l/2 + ) = 0.95 ±0.10 GeV, (14) 
for the same doublet when the currents @, @ with the derivative is used. 

A(3/2 + ) = 0.82 ±0.10 GeV, (15) 

for the doublet (l + ,2 + ). These results shall be used in the following sections. 
In the following sections we also need the values of f 's. 

/ +3/2 = 0.19 ± 0.03 GeV 5/2 , (16) 

f + l/2 = 0.37 ± 0.06 GeV 5/2 , (17) 

/ +1/2 = -(0.40 ± 0.06) GeV 3/2 . (18) 

Note it is impossible to determine the absolute sign of f 's from QCD sum rule approach. 
From the equation of motion we can derive Qj: 

f+,1/2 = -A + ,|/+,i/2 ■ (19) 
In this work we adopt the convention of the positive f' + 1//2 , /+,3/2- 
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3 The sum rules at the 0(1/ mq) order 



Inserting the heavy meson eigen-state of the Hamiltonian up to the order O{l/mo) into 
the correlator 



'7 d 4 x^(0\T^^(x)j]^(0)) |0> , 

the pole term on the hadron side becomes 

(f + 5ff f 25mp 2f5f 



n(w) 



pole 



2(A + 5m) - u 2K-uu (2k - lu) 2 2A - lu 



(20) 



(21) 



where 5m and 5f are of the order 0(1/ mo). 

To extract 5m in (|2lD we follow the approach of and consider the three-point corre- 
lation functions 



jf ||":-".- , , ,' 



<i\n0,i<»- ;A '- ? '(ov (./;',:;" i.d o(o) jJJv- ' (//)) |o> , (22) 



where O = K or 5 is the kinetic energy or the spin dependent term of the 0(1/ mq) 
lagrangian of the HQET. The scalar function corresponding to (|22]) can be represented as 
the double dispersion integral 



/ 1 r 
5 U(uj,u}') = — J 

The pole parts for them are 



p (s, s')dsds' 
(s — oj)(s' — uj') 



(23) 



5 K U(u,LU ) po le 

5 s ll(uj,uj)p i e 



PK 



(2A-u;)(2A-a/) 
d M P^ 



7T + 



pG K (u') , pG K (u) 



(2A-w)(2A-u/) 



77- + d M f 



2k -lu 2A-u/ ' 
G 5 (u/) , G 5 (u;) 



2A - 2A - w' 



where 



K j,P,h = (j, p ,3i\hv (iD ± ) 2 h v \j,P,ji) , 

2d M T, jjPijl = (j, P, ji\h v ga^G^Klj, P, j t ) , 

dm = d jjjv djt-ifijt = 2ji + 2, dj l+ i/2^ = —2ji. 

Let lu = lu' in (^4j) and fl25|), and compare it with fl2~T|) one obtains|| 



<5m 



4mr 



■(A + dMC mag Y,) . 



(24) 
(25) 



(26) 
(27) 
(28) 



(29) 



The simple pole term in ( p4"|) and (|25| ) comes from the region in which s(s') = 2A and s(s') 
is at the pole for a radial excited state or in the continuum. These terms are suppressed by 
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making double Borel transformation for bo th ui and u'. One obtains thus the sum rules 
for K and S as 

f Ke -^/T = r r dwdw'e-^'y^pAu^') , (30) 

Jo Jo 

/*£ e - 2K ' T = r r dujdu'e-^ + ^l 2T Ps {uj, u') , (31) 
Jo Jo 

where the spectral densities are obtained from straightforward calculations in HQET. 

Confining us to the leading order of perturbation and the operators with dimension 

1 + 

D < 5 in OPE, we find for the jf = — doublet, 

f K e -2A/T = 3 r ^ e-«' T du + JL (a s GG) T 3 , (32) 
2'ti z Jo 2rix 

/ 2 S e" 2A/T = r r P(s, s') e~^dsds' + — (a s GG) T 3 , (33) 

Jo Jo 487T 

with p(s, s') = -^ as( ?P C ma gSs' {s' 2 {3s — s')9{s — s') + (s <-> s')} when the currents © and 
(IBD are used and 



f Ke -2A/T = * f u < e -./T du _ _L {asGG) T- 6 -ml (qq) , (34) 

2°7T Z JO 2 7T O 

/ 2 S e~ 2A/T = jf* 5') e-^rf^' + ^{oc s GG) r+imj (g ff > (35) 



with p(s,s') = ^ "^3 C'mag{g /2 (3g — s')0(s — s') + (s <-> s')} when the currents (0) and 

3 + 

([|) with no extra derivative are used. For jf = — doublet, a straightforward calculation 
yields 



fK^lr = __L f J ^ + ^(a.GG) T* , (36) 
f 2 E e -2A/r = / / , a e -4rdsds' + -=-(a s GG) T 3 . (37) 

Jo Jo 727T 



with p(s, s') = ^ ^pC mag {s'\s - ls')9(s - s') + (s^ s')}. Combining fD-fZl) with 
(0p-(O) we can obtain sum rules for K and £ in the three cases. 

The spin-symmetry violating term S not only causes splitting of masses in the same 
doublet, but also causes mixing of states with the same j, P but different ji. This mixing 
is characterized by the matrix element 

(3,P,J + \\ K g -°»»G^h v \j,P,j - i) = -2 m + 4J p ) . (38) 

This quantity can be extracted from the correlator 

i 2 J d 4 xd%e tk -*- M -y(0\T (j hP , + ±{x) h v 9 -a, u G^h v (0) J.^ j//)) |0> , (39) 
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the double pole term of which is 



With the methods similar to that used above one can obtain the sum rule of m + _(J p ). 
For two 1 + states we find 

rm a(l + ) f' + ±Uie~ R+h/T e~ K+ * /T = jT jT P (s, s') e~^dsds' + ^ (a s GG) T 3 .(41) 

with p(s, s') = ^^¥ 1 C mag {ss li 6(s - s') + s 3 (6s' 2 - 8ss' + 3s 2 )9(s' - s)} for the currents 
(D and ©. 

The sum rules (32)- (37) and (41) were obtained in [|J except that all loop corrections 
of the order a s were neglected there. In this work we have included the perturbative terms 
of the order a s in the sum rules for £ and ^1/2,3/2 which are leading perturbative terms 
for these quantities. 

The masses of two 1 + states are obtained by diagonalizing the mass matrix 



^A, I - 7^— {K, 1 - C mag E + 1) — C mag mi a(l + ) ^ 

+ '2 4m n + '2 y + '2 y 4m n 2 ' 2 



1 1 

^ C maff m fj i (1+) A +if - — (K + , + 5 C mag 



(42) 



where C ma9 = , /3o — H — 111 our analysis we keep four active flavors with 

the QCD parameter Aqcd = 220MeV. 

We can minimize the dependence of the sum rules obtained in Section |3| on /, A and 
u c through dividing the sum rules in (|3~2D-(|3~7D by the sum rules in (^)-([TT|) respectively. 
The seemingly stability at large T values is not useful, since in this region the sum rules is 
strongly contaminated by higher resonance states. This continuum model contamination 
problem is quite severe also in the case of sum rules analysis of ground state heavy meson 
to 1/rriQ H [n§. It originates from the high power dependence of the spectral densities. 
Moreover the stability plateau for the three-point sum rules does not necessarily coincide 
with that of two-point sum rules. We follow and use the same working windows as those 
of two-point sum rules in the evaluation of 0(l/rriQ) corrections. 

Taking the working regions of the sum rules as that for the two-piont functions for the 
three cases, we obtain the set of values for K as following: 

K = -2.30 ± 0.40 GeV 2 , (43) 
for the doublet (0 + , 1 + ) when the currents (|3|), (Q) without the derivative are used, 

K = -1.70 ±0.40 GeV 2 , (44) 
for the same doublet when the currents (j^), ((^) with the derivative are used, and 

K = -1.60 ±0.40 GeV 2 , (45) 
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for the doublet (l + ,2+). 

The results for E for the doublet (0 + , 1 + ) in the working regions of the sum rules are 

E = 0.22 ± 0.04 GeV 2 , (46) 

when the currents (|j) without the derivative are used and 

E = 0.18 ±0.03 GeV 2 , (47) 

when the currents ©, @ with the derivative are used. For the doublet (1 + , 2 + ), the value 
of E is found to be 

E = 0.03 ± 0.01 GeV 2 . (48) 

Dividing the sum rule in ( |4T| ) by the sum rules in ( |T0| ) and (pTTj) , we also obtain the 
expression for mi i(l + ). Using the working windows as u c = 2.8 — 3.2 GeV and T = 
0.53 — 0.73 GeV, we obtain 



mi |(l + ) = 0.08 ±0.02 GeV 2 . (49) 



The errors quoted above reflect only the variations with the Borel parameter T and the 
continuum threshold uj c within the working windows. They do not include other intrinsic 
errors in QCD sum rules in HQET. The corrections to the masses due to the mixing 
are formally of the order —r. By diagonalizing the mass matrix we found that they are 

m Q 

numerically negligible. 

Since the experimental data is not yet available for the excited B mesons, we present 
our results for the D system assuming the HQET is good enough for the excited D meson 
too. For the doublet (0 + , 1 + ), we have 

- {m D * + 3m Di ) = m c + A + — [(0.57 ± 0.10) GeV 2 ] , (50) 

m D > -m D * = — [(0.22 ± 0.04) GeV 2 ] , (51) 
1 m c 

when the currents @ and @ are used and 

- (m D * + 3m D , ) = m c + A + — [(0.42 ± 0.10) GeV 2 ] , (52) 
4 1 m c 

mjy -m D * = — [(0.18 ± 0.03) GeV 2 ] , (53) 

1 >>i 



Til. 



c 



when the currents (||) and @ are used. As for the doublet (l + ,2 + ), the result is 

- {3m Dl + 5m fl .) = m c + A + — [(0.40 ± 0.10) GeV 2 ] , (54) 
8 2 m c 

m D * -m Dl = — [(0.06 ± 0.02) GeV 2 ] . (55) 

2 m c 

The renormalization coefficient C mag of the chromomagnetic operator for charmed meson 
in the above formulas is reasonably neglected. The results for B system are obtained by 



8 



replacing m c by m b and multiplying the right hand sides of (51), (53) and (55) by 0.8 since 
Cmag ~ 0.8 for B system. 

The results in ( [54]) with the A value in (|15D and m c = 1.4GeV agree reasonably with 
experimental data, though the uncertainties, especially for A, are somewhat large. 

In the following sections we will need the values of r x , 3 = 5f ± , 3 / f x + 3 , r x _ 1 = 
5f 1 _ 1 //j _ 1 and 5/i 3 , where 5/ is the correction to / due to the K, and S terms in the 
Lagrangian at the order of 0(1/ tuq) and Sfi 3 is defined as: 

(0|J 1)+! x|l,+,^> = Sf u r] lj+)l (56) 

which is also of the order I/tuq. 

Omitting the indices j, P,ji, r can be written in the form: 

Gk . , Gy, . s 

' = 2^ + *<2^ ' (57) 

where Gk and G-% are defined as: 

(0\i / d*xK.(x)J jiPjl (0)\j,p, ji) = ffojiGicVjipj, , (58) 



(0|i j d xS(x)J jiPijl (0)\j,p,ji) = fj, P , jl d M G Tl r] j ^ jl . (59) 
The mixing coupling constant 8 fx 3 = Gi zjlmQ where G\ 3 is defined as: 

/ ^5(x)J lf+ , J (0)|l > + ) 2) = Gj ii%|+ii . (60) 

r x _ j 1 for the ground state 1~ has been calculated in |19|, ^0|. Notice that bj x _\ calcu- 
lated in 0, |20| contains an additional term proportional to A/tuq. This term comes from 
the 0(1 /mo) correction in the expansion of the physical vector current Q'j^q responsible 
to the leptonic decay of B*, where Q is the field in full QCD. For our purpose instead the 
matrix elements of the interpolating currents Jj, p j t defined in ([!]) and @ at the leading 
order are needed. Therefore this term is not included in fl5"T|). ^From the results of pOf , 
we know i = —(0.82 ± 0.30)/mg for 1~ and r _ 1 = —(0.74 ± 0.30)/mg for 0~ heavy 
meson at the leading order of a s [pOR . 



We follow the method used in [201 to derive the Gk, Gt, and Gi 3 for the B\ meson 

L| — 1 ' 2 ' 2 

and B\ doublet. Letting k = k' and u = oj' in (p2])-( 25|) and (|39|)-(|40D we arrive at the 
following sum rules: 

fl*e-^{2G K + £) = --jL r 8 *e-*d8 + l.p.(?)T*, (61) 
^'2 1 64vt z Jo 96 71 

„9 iA +'l , „ S. a s r u <= r.25 4 , s, » , 1 ,at a _ 9 . . . 

^—(2^ + -) = ^ s =(_-_ 1 n-)e- Tds + -<^G^, (62) 
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/+ 3/i ie * (-H^ + -^±-) = i— f / s 5 (- - - In -)e~Tds + ^-(^G 2 )T 2 



' 2 



f + i T 32tt 3 7o v 81 9 72 N vr 

(63) 

where we take /i = 2T as in the evaluation of K and S. 

In order to obtain Gk, Gy, and Gi | we subtract (0), ( ]62| ) and (|63"D by the expressions 
of (36), (37) and (41) after they are divided by T. The dependence of Gk, G% and 
Gi a//, i on the parameter T and c<j c is shown in Fig. 1-3 with uo c = 3.2,3.0,2.8 GeV. 
In the numerical analysis we choose the same stability region of (T, lu c ) as those from the 
mass sum rules. In this region there exists stability plateu for the sum rules (|B"ID-(|6"2D. For 
the sum rule ( |63[ ) the gluon condensate contribution is much bigger than the perturbative 
one, which results in the curves in Fig. 3. 

Numerically we have: G K = -(1.0±0.20±0.25)GeV, G E = (0.013±0.005±0.002)GeV, 
Here the first error comes from the variation with u c and the second error is due to the 
variation with T in the working region. Giz/f + i = (0.023 ± 0.003 ± 0.010)GeV. Finally 
we arrive at r 1+ 3 = -(0.52±0.22)/m Q , r 2 ' + % = -(0.57±0.25)/m Q , (J/i >2 = (4.3 ±2.0) x 
10- 3 /m Q GeV 5 / 2 . 2 



4 Mixing of the Two 1+ States 



Diagonalizing the mixing mass matrix ([42]) we arrive at the physical eigenstates for the 
two 1 + states up to the first order of I/tuq. 

|Si>=cos0|l,+,i>-sin0|l J + J |) J (64) 

|fi 1 ) = sin^|l,+,i)+cos^|l,+J), (65) 
where 9 is the mixing angle, which is determined by the folowing equation: 

tciii 2$ = 2 ' 2 f 66) 

4m Q (A +i i - A +)f ) + (K + i - K + s) + C mag (i: + ^ + 5E +if ) ' 

Numerically we get 9 = —(0.10 ± 0.05) for the D system with m c = 1.4GeV and 9 = 
-(0.03 ± 0.015) for the B system with m b = 4.8GeV. 

The corresponding interpolating currents which annihiliate one of these two states and 
do not couple to the other up to the first order of 1/rriQ can be written as 

■' fli = V,*+« J i,+,§' (67) 

J Bl = J h+> 3 + PJ 1>+ i ■ (68) 

In the discussion of the mixing of the two 1 + states we always use the current (^) with 
derivatives. The mixing parameters a, ft are determined from the orthonormal conditions 
of the physical 1 + states: 

<0|J^|Si) = 0, (69) 



10 



(0\J Bl \B[) = 0. (70) 
Solving the above equations we obtain 

f , , i tan 9 + 5 f i 3 
a = ^ (71) 

f i , 3 tan# — (5 f 3 i 
/? = (72) 

With the values of f's in flU)-®, <Vi,a and tan # we S et « = °- 22 ± - 10 > Z 5 = 
-(0.06±0.03) for D[, D l and a = 0.065 ±0.0 2 32, (3 = -(0.02±0.01) for B[, B l respectively. 



5 Pionic decay amplitudes of B\ up to the order of 

0(l/m Q ) 

In this section we shall use the light-cone QCD sum rules (LCQSR) [§] to calculate the 
partial widths of the pionic decay processes of B\ and B\. The decay width of (l + ,2 + ), 
(0 + , 1 + ) doublets have been calculated at the leading order in HQET with conventional 
QCD sum rules |J. Different from the conventional QCD sum rules, which is based on the 
short distance operator product expansion (OPE), LCQSR is based on the OPE on the 
light cone which is the expansion over the twists of the operators. 

Denote the doublet (l + ,2 + ) by {B\^B^). Let us consider first the decay process B\ — > 
B*7i. ^From (^) the decay amplitude can be written as 

M(B 1 -> B*ir) = cos 9M{- ) + sin#M(- ) , (73) 

where M(| + ) and M(| + ) denote contributions from the states ji = | + an d 3i = \ + 
respectively. Up to the order of (9(l/mg) we can confine us to the leading order for M(i + ) 

and up to the next-to-leading order for for M(| + ) in the l/m<g expansion. 

^From covariance and conservation of the angular momentum of the light component 
in the rriQ — > oo limit, the leading order amplitudes have the general form: 

M (B 1 -> B*tt) = M (^ + ) = / e; V M<lt ~ \gF<g)g{Bx, B*) , (74) 

M^ + ) = Ie*.r ] g'(B' 1 ,B*), (75) 

corresponding to the D-wave and S-wave for the pion respectively. Here r/^ and are 
polarization vectors for the states 1 + and 1~ respectively. q tfJ , — g M — v ■ qv^. I = v^2, 1 for 
charged and neutral pion respectively. 

The 0{l/mQ) correction Mi(| + ) arises from the the terms /C and S in the Lagrangian 
in HQET. The operator K, contributes only to the D-wave amplitude, while S contributes 
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to both the D-wave and S-wave amplitudes. Therefore the total 0(l/mg) correction of 
the decay amplitude M(Bi —>■ B*ir) has the general form: 

M l (B 1 - B*tt) = M u + Mis — I e;vA{Q?<lt ~ ^f<7t + stgxs} ■ (76) 
In order to get a consistent expansion over 1 / rriQ we consider the correlator 

/ d\ e-^(n(q)\T (j^(0)4>)) |0) , (77) 

where Jb 1 is the interpolating current ( |S8"D for B\. It is implied that the 0{l/mq) terms 
in the Lagrangian of HQET are included in the calculation of (|77|). 

The pole term of (|77|) can be expressed up to the order 0(1/itiq) as: 

+ Sf_i)(f + 3 + 6f + *)(M + Mu + MuT*_ 
[2(A_ i + 8m_x) - cj'][2(A +j 3 + 5m + z) - to] 



2(A +) | + 5m +) |) - io 2(A_ j i + 5m_ i ) - u' 
which can be expanded as: 

f-,lU,l M f , /-,s/ + ,§(ri,- f4 +r 1>+ ,|)^ + /_,|/ +1 |(M W + M ls )^ 



(2A_i -u/)(2A +j | -w) (2A_i -o/)(2A +) s - w) (2A_i - u/)(2A +> 3 - oo 

/_ i/ + zlbmxMf f_ if, s25m3Mf 



2 



+ ■■■ ,(79) 



(2A_ ) i - w') 2 (2A +j | - u) (2A_ ) i - cu')(2A +j | - cj) 5 

where M t - , z = 0, 1, is defined by Mj = ej^M"' 3 , and r r _ i, r 1+ s is defined in Section 
(§• 

On the other hand ([77]) can be expanded as: 

J d 4 xe- tt -<7r(g)|T(jj;_ j i(0)^ +ti (x)) |0> 
+/? / rf 4 x e-* fc -<7r(g)|r(^_ i i(0)^. ti (x)) |0> 

| d 4 xA e^V^^T^^^ |0) 

+z | d 4 xd 4 y e- ik - x e lk '-y(7c(q)\T (j^ ^^h^DtfKiO) j£ f (x)^ |0> , (80) 

where the first term is the correlator at the leading order of I/ttlq expansion, the second 
term is due to the mixing of two 1 + currents, the third and fourth term arises from the 
C(l/mg) corrections to the transition j e = | + to ji = | due to the O{l/mo) terms of 
the Lagrangian in the order of Oil/mg). These correlators will be calculated with OPE 
and the leading order Lagrangian of HQET. 
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6 The sum rules for the coupling constants g, g' etc 

We consider the correlators appearing in ([BOf). They have the general form: 

J d 4 x e-^(n(q)\T (^_,i(0) (x)) |0> = (gft? - ^f\ 2 ) G(u;,a/) , (81) 

/ rf 4 x e-^(7r(g)|T (j;„ ,(0)j;^(x)) |0) = tf*G>(u,<J) , (82) 

• y rf 4 iA e- ifc -e^(7r( g )|T (j^G/), jK(0)a ■ Gh v (0), ^(z)) |0> (83) 



^r\ 2 )Gf d (^^)+^Gf s (^,^)) , 



z J d 4 xd 4 y e- lk -*e lk '-y(n(q)\T ^ |0> 

= ^ - ^f\ 2 ) a/) , (84) 

where k! = k — q, u = 2v ■ k, uj = 2v ■ k' and q 2 = m 2 « 0. The suffices K, and 5 denote 
the contributions from the operators from /C and S respectively. The forms of the right 
hand side of fl8~I])- (jS2l) are determined by that a and (3 are transverse indices, x — y = vt 
on the heavy quark propagator and the conservation of angular momentum of the light 
component. 

For deriving QCD sum rules we first calculate the correlator (|8T|) , which is equal to 

Y C dt I ^^(-v-^W + tMDa - |^)Mff)l«(0)d(x)|0>} , (85) 
By the operator expansion on the light-cone the matrix element of the nonlocal opera- 



tors between the vacuum and pion state in (85) defines the two particle wave function of 



the pion. Up to twist four the Dirac components of this wave function can be written as 



< n(q)\d(x) W u(0)\0 > = ~^U% / du e luqx {^{u) + x 2 9l (u) + 0(x 4 )) 

J 

x 2 q^, f 1 



+ / due m « x g 2 (u), (86) 

qx Jo 

- f vn 2 

< n(q)\d(x)i^u(0)\0 > = Jw w / du e mqx (p P {u) , (87) 



m u + m d Jo 

— f in 2 

< 7r(g)|d(a;)<V7 5 w(0)|0 > = i{q„x v - q u xj . J \* / du e^cp^u) . (88) 

6{m u + m d ) Jo 

The wave function (p n is associated with the leading twist two operator, g\ and g 2 corre- 
spond to twist four operators, and tpp and (p a to twist three ones. Due to the choice of 
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the gauge x^A^(x) = 0, the path-ordered gauge factor Pexp (ig s Jq dux^A^ux)) has been 
omitted. The coefficient in front of the r.h.s. of eqs. (|8T|), (|88| ) can be written in terms of the 

ml _2_ 

m u + m d fl 



light quark condensate < uu > using the PCAC relation: ^ = - — = — — < uu>. 

' J- z 



From (p5[ ), (j86p and (|88|) we arrive at 



G(co,u') = / dt / due^-e^-uW 7T (u)+t 2 g 1 (u) + g 2 (u)+ - fjL% < Pa ( u )}+- ■ ■ . 

8 jo jo q ■ v o 

(89) 

For large euclidean values of uj and u' this integral is dominated by the region of small t, 
therefore it can be approximated by the first a few terms. 
Similarly we have: 

G>, to') = -{qlU / °° dt Jo 1 due^f e^u{^{u) + t 2 9l (u) + ^g 2 (u) + § fi^u)} + ■■■ 
= +±(co - lu') 2 UJ™ dtj, 1 due^~ u ^e iu ^u{^(u) + t 2 g x {u) + f v g 2 {u) + ^<p a (u)} + ■ -(-99) 

with the current (|6]) for ji = | + state. Or 

j poo rl , i, 

G'{u,J) = -U / dt / due^-^-e^-{^ P {n)-{q-v)[^{u)+t 2 gi {u)}} + --- . (91) 
4 Jo Jo 

with the current (||). 

After double Borel transformation with the variables uj and uj' the single-pole terms 
in (|78| ) are eliminated. Using (0)) (0)> (H) an d subtracting the continuum contribution 
which is modeled by the dispersion integral in region u,u' > u c , we arrive at: 

y/E A -,j +A +4 _^ 4 4 1 

= T {tt ^(tt )T(l-e T )-^Wo#i(wo) + yG 2 (w ) + 3/VWoy^(wo)} , 

(92) 

with = ^ = 92MeV, tt = T ^ T2 , T = j^pjr , where Ti, T 2 are the Borel parameters and 
G 2 («o) = Io° u 92(u)du. In obtaining ([52|) we have used the Borel transformation formula: 
BZe°» = 6(a- ±). 

Similarly, for the coupling constant we have: 

= R e ~'" T (w^W^-Mf ) - 4^!( U )T/ (f ) + |^ CT (n)T 2 /i(f )) |u=« 

n k 

with the current (|6]), where f n (x) = 1 — e~ x fr is the factor used to subtract the con- 

k=0 

tinuum. The derivative in ( |93D arises from the factor q 2 in (0). We have used integration 
by parts to absorb the factors q ■ v and (q ■ v) 2 . In this way we arrive at the simple form 
after double Borel transformation. 
With the current ffl) we have 



9'f-if + i = - A F*e -WWTVilS + 2^ P ( Mo )T/ (^) + ^(t*,)} , (94) 



A_ 1+ A' 

14 



where ip' w (u ), g[(u ) are the first derivatives of <p n (u), gi(u) at u = u . 

For the 0(1/ mo) correction to the decay amplitude for (j e — |) — > (jg — |) + 7r, we 
need to consider the correlators (|83]) and (plf ). In calculations of these correlators enter 
the pion matrix elements of quark-gluon operators, which are parameterized in terms of 
pion wave functions defined in ||: 

< n(q)\d(x)a a p'j 5 g s G^(ux)u(0)\0 >= 

if^q^qaQup - q v q a 9tii3) - (q^qpgua - q^qpg^)] / Von ^ 3n (ai)e iqx{ai+va3) , (95) 



< n(q)\d(x)j^ 5 g s G a p(vx)u(0)\0 >= 

q ■ x 



f*[qti(9*n- ^rz) -?«(. 



q ■ x 



gpn 



Vaiip ± (ai)e 



iqx(oL\-\-voi.z) 



Q 

+U— —(q a X/3 - qpx a ) I Va^Aa^e 
q ■ x J 

< Ti(q)\d(x)-i^g s G af3 (vx)u(Q)Y} >- 

•^a.q^, 



and 



a iqx{a\+vai) 



(96) 



if* 



qp(ga» 



q ■ x 



q a [gpn 



q ■ x 



+ifn— l! -(qaX(3 ~ q/3Xa) / X>Oi£|| (c^e 



q ■ x 



Vaiip ± (ai)e iqx{ai+va3) 

iqx(cti+vct3) 



(97) 



The operator G a p is the dual of G a /3- G a p = \e a ^ p G &p \ Van is defined as Vaii = 
da>ida>2da>38(l — a± — a 2 — 0:3). The function (p 37r is of twist three, while all the wave 
functions appearing in eqs.(|9~6"|), (|97l) are of twist four. The wave functions <p(xi, /x) (/x is 
the renormalization point) describe the distribution in longitudinal momenta inside the 
pion, the parameters x« (J2i — 1) representing the fractions of the longitudinal momen- 
tum carried by the quark and the antiquark. 

The wave function normalizations immediately follow from the definitions d86|)-(p7|): 
J^du (p n (u) = J^du <p a (u) = 1, fgdu gi(u) = S 2 /12, J Va^^ai) = f Va^^ai) = 0, 
/ Vcti(p±(ai) = — J Va.i{p\\(oii) = S 2 /3, with the parameter 5 defined by the matrix element: 
< i^(q)\dg s G a ^ a u\Q >= w^/vr^- 

For sum rules for the 1/tuq corrections we shall compare the double Borel transform of 
(|7|) with those of ( p3| ) and (|4j) for the D-wave and S-wave part seperately. It is convenient 
to confine to be longitudinal, i.e., k l = 0. For the D-wave amplitude we have: 



VE 1 

12 rriQ 



t-^ttt / ' dt x dt 2 Va % oi 2 \-i^(q- v) V ^(a/) + F^ L (a/) + ^^(a/)} 

e ikvt e -iqv[a3t2+a2(ti+t2)} 



G?>,u/) = |Gf>,u/), 



G«(u, u>) = f ^ f °° dhdh Jo 1 due^r-^ e 

+(ti + h) 2 gi(u 



2 (k - uq) 2 u{ip n (u) 



^lg 2 ( U ) + ^±^^ a (u)} + 



q-v 



(98) 
(99) 

(100) 
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where (k — uq) 2 = —^(uj — uj') 2 



4 

After double Borel transformation we have: 



8m i +<5m 3 

[(n _ : i + r l+ z)g + g ld V— -9\f- 

A_ 1 +A 3 



+f 5j / ^iat T /o(f ) + T)} , (101) 

where the function K(uq,T) arises from the kinetic operator /C, which is defined as: 



K(u ,T) = ^ («V,(«)T 2 /i(f)-4« 3 9iW + |^V.HT/ (f )) \u- 



+A£(u 3 g 2 (u))\ u=U0 . (102) 

In order to extract g\& we have to use the operator -§f\T x ( |101| )}. 
For the S-wave part we find 

I (a^M)) Tfo(f) 



+I F A£ (u^(u)T 3 f 2 (f) - 4u gi (u)Tf (f) + \^u Va {u)T 2 h\ 



OJc > 

T J J \u=uq 



+4±(ug 2 (u))\ u=U0 Tf (f)}, (103) 



where f3 is defined in section f| and its value is found through (|72]) 



7 Sum rules for the pionic decay amplitudes of B\ up 
to the order of 0(l/mg) 

In the following we shall consider the 0{l/rriQ) correction to the pionic decay amplitude 
of B\ meson. There are two decay channels. The decay amplitudes are: 

M(B* — ► Bit) — I V^qt[9 a (B;, B) + 9 <{{B*, B)] , (104) 
M(B* 2 - B*tt) = J ie aPtjp e* a v^q p t q t ,[g\Bl, B*) + ^(5*, B*)} , (105) 

where r/^ is the polarization tensor of B\ meson, g a , g b are the coupling constants at the 
leading order of l/mq and gf, g\ is the Oil/mq) correction to g a , g b . It can be shown 
fl7|l by combining heavy quark symmetry and chiral symmetry that at the leading order 
of I/nig the coupling constants in ([73]), ( |104j ) and ( |105| ) satisfy 

g a (B*,B) = g b (B*,B*) = ^g(B 1 ,B*) . (106) 

Angular momentum conservation requires that both operators K, and S contribute to the 
D-wave amplitude only, which is confirmed through explicit calculation. Moreover we do 
not consider mixing effects for B\ meson, which is expected to be small. 
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Since the detailed calculation is very similar to that for Bi meson, we present the final 
results only: 

8m i +8m 3 

[(r 0| _ fi + r 2)+;| )^ + g$ - -i^^f^f^ 

= e "K +4 {^-JVa^ ± ( ai ) - 1^1,(0,)] 
-4 fe/^ii («2^K))T/ (f ) + g%,T)} , (107) 



Jtn 1 +om 3 

+ r 2, + ,§)^ + g\ - -V-Vl/-,*/*,* 

= — {fVai^niai) + \K{u Q ,T)} , (108) 

where the function K(uq,T) is defined in ( |102| ). 

Note that the corrections to the leading order coupling constants from the kinetic 
operator K satisfy the same relation as ( |106|) , which is an expected result since K, does not 
depend on the spin. 



8 Numerical analysis of the sum rules 

In the numerical analysis we shall take 7\ = T 2 = 2T and Uq = 1/2 as done in || and ||. 
Although the initial and final heavy meson states are different in their masses, this seems 
to be a reasonable approximation in view of the large value of T\ and T 2 used below. 

Note that our sum rules depend on the pion wave functions and the integrals or deriva- 
tives of them at the point uq = |. Since 5 2 is numerically small, the uncertanty due to 
the integral term G 2 (ito) is insignificant. The values of the wave functions at Uq have been 
determined quite well and their dependence on the renormalization scale is very weak M. 
The first derivatives at uq are zero. The second derivatives are very sensitive to the detailed 
shape of the wave functions. Usually only the lowest few moments of the wave functions 
can be determined by QCD sum rules. There are many wave functions satisfying these 
constraints from moments. With the form of wave functions in ||, the final sum rules are 
sensitive to the renormalization scale and the coefficients cij in ip n (u) due to its large value. 
This fact indicates that the not- well-known nonperturbative contribution from condensates 
in the QCD vacuum and the even higher twist wave functions might be important. In order 
to estimate the corrections from the above sources we present the detailed expressions of 
the pion wave functions and their second derivatives in order to show the dependence on 
the parameters u>ij in them. 

There are many discussions about <p n (u) in literature [0, [13], |[ [14|, [15], [16| . The model 



wave function for <p w (u) based on the QCD sum rule approach was given in as: 

(p v (u,n) = 6uu (l + a 2 (/i)-[5(w - u) 2 - 1] + a 4 (/x) — [21 (it - uf - lA(u - u) 2 + 1]) . 

V 2 8 / 

(109) 

Yet the authors pointed out that the oscillations of ( |109|) around uq = \ is unphysical, which 
is due to the truncation of the series and keeping only the first a few terms when <p n (u) is 
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expanded over Gegenbauer polynominals. In [15] it was stressed that: (1) the expansion 
over Gegenbauer polynominals converges very slowly as can be seen from the large value 
of a 2 and a 4 ; (2) any oscillating wave function is not physical, since no detached scale is 



seen to govern such oscillations. Recently Mikhailov and Radyushkin [I6| reanalysed the 
QCD sum rules for <p n (u) taking into account the nonlocality of the condensates. They 
suggested the following wave function: 



<p„(u) = -Ju(l-u) , (110) 

7T v 



which is close to the asymptotic form and the suggested pion wave function in ]T5 

(771 2 \ 
~ %(3 2 u(l - u) ) + W + trf ~ 2 ) u ( l ~ u )} > ( ln ) 

where N = 4.53, p. = ||, p, = ± (\m v + |m p ) = || with m = 330MeV, rh = 620MeV, and 
(3 = 320MeV. The Brodsky-Huang-Lepage wave function is also close to the asymptotic 
form except for the different behavior at end points Jl~| : 

/ M 2 \ 

cp^u) = 9.05^(1 - u)exp -^-p. r , (112) 

\ Sa z u(l — u) J 

where M = 289MeV and a = 385MeV. The asymptotic form of ip n (u) reads: 

^ m {u) = 6w(l - u) . (113) 

For the other pion wave functions we use the results given in ||. The detailed expres- 
sions of the pion wave functions relevant in our calculation are: 

<P3n((*i) = 360aia 2 a 2 [l + ^i,o^(7a 3 - 3) + ^ 2>0 (2 - 4aia 2 - 8a 3 + 8a 2 ) 

+wi,i(3aiot 2 - 2a 3 + 3a|)] , (114) 

<p P {u,n) = I + B 2 {p)\[3{u - u) 2 - 1] + 5 4 (^[35(w-w) 4 -30(w-w) 2 + 3], (115) 

2 8 

(p a (u,fi) = 6uu (^l + C 2 (/i)^[5(n-n) 2 - l] + C 4 (/i)y[21(n-M) 4 -14(n-n) 2 + l]^ , 

(116) 

<p ± (on) = 30S 2 ( ai - a 2 )a 2 [i + 2e(l - 2« 3 )] , (117) 
011 (a,) = -1205 2 ai a 2 a 3 [^ + e(l - 2a 3 )] , (118) 

5 1 6 

gi (u) = -5 2 u 2 u 2 + -e5 2 [uu(2 + 13uu) + 10M 3 lnw(2 - 3m + -u 2 ) 
2 2 5 

+I0u 3 lnu(2-3u + -u 2 ) , (119) 
5 
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10 

g 2 (u) = —5 2 uu{u - u) , (120) 

G a («) = -f5 a « 3 (|« a -|«+|) ) (121) 

where u = 1 — u, B 2 = 30R, B 4 = §i?(4u; 2i o — ui t i — 2^0), C 2 = R(5 — 5^1,0) > C4 — 

^i?(4tt; 2i o — u>x,i) and i? = The coefficients f^-, Wi,k 01) $ 2 @ an< l e have been 
determined from QCD sum rules. 

We shall take the scale /1 = %uj c = 2.4GeV, at which the various parameters are: 

a 2 = 0.35, a 4 = 0.18, u lfi = -2.18, u 2fi = 8.12, u 1A = -2.59, f 3n = 0.0026GeV 2 , 

= 2.02GeV, S 2 = 0.17GeV 2 , e = 0.36. 

3 3 15 

= -{I- -a 2 + —ai), (122) 

M«o) = 1-~B 2 + ^ 4 , (123) 

^(«o) = §(l-|c7 2 + yC4), (124) 

3iW = ( 4 - 0.037e)5 2 , (125) 



32 

~288' 



G 2 (u ) = - (126) 



225 

<p'l(uo) = -12(1 - 9a 2 + — a 4 ), (127) 

o 
22 F) 

^(uo) = -12(1 - 9C 2 + — C 4 ), (128) 

o 

<tfM = -(2.5 + 7.4e)5 2 , (129) 
fl£(«o) = ~<* 2 , (130) 

/p a .% ± ( a .) = _^L(l + 2e ), (131) 
j a 3 12 

Da,%ll(a l ) = -(^ + 6)5 2 , (132) 
«3 3 

1 d r / m 4 8 . . 

Va i 3 [Ot2<P3A a i)\ = 6 - 2wi o - -^2,0 + , (133) 

a3 aa 2 7 7 

where the parameters a.j, Cj, tjjj and e indicate the deviation from the asymtotic form 
of pion wave functions. 

At Mo = | the values of the functions appearing in (|9l^)-( |T0~T| ) are: <p n (uo) = 1.22, 
1.273, 1.25, 1.71, 1.5 for the wave functions (|109|) - (|113|) resepctively, <pp(uo) = 1.07, 
ip a (u ) = 1.55, gi{u Q ) = 0.022GeV 2 , g 2 (u ) = and G 2 (u ) = -4.1 x 10" 3 GeV 2 . Their first 
derivatives satisfy: (p' n (u ) = ip' P (uo) = ^(uq) = g[(u ) = G' 2 (u ) = 0. And their second 
derivatives are: (p"(u ) = -35, -5.09, -0.02, -17.50, -12 for (PTJP|)- (|TT5|) resepctively, 
ipZ(u ) = -17, g'{{u ) = -0.88GeV 2 , g' 2 (u ) = 0.28GeV. 



/ 



19 



Note the theoretical side of the sum rules ( |101| ), (|107|) and (|108|) depends on f w (u) only 
through the combination [u 3 <£v(u)]" o = 3ip n (uo) + ^ip"(u ) in K{u Q) T). Its value is 3.18, 
3.75, 2.93 and 3 for ( |110| )- (|113|) resepctively, which is relatively close to each other. We 
shall present numerical results using the wave function ( |1 1 2| ) . 

The dependence on the Borel parameter T of f_ if + ig, f_ if + ig', mQ/_i/ + 3j ls , 

m Q f_if + 3 ((r lt _i+r lj+ *)g + g ld ), m Q f_if + 3 ((r 0> _ i + r 2+ *)g a + gf), and 

"^qZ-^Z+j i + r 2+ 3)g b + gfj are shown in Fig. 4-10 with u c = 3.2,3.0,2.8 GeV 

and the pion wave function ( |112]) . Requiring that the higher twist contribution is less than 
30% of the whole sum rule, we get the lower limit for the Borel parameter T, T > l.OGeV. 
Requiring that the continuum contribution is less than 40% of the whole sum rule, we get 
the upper limit, T < 2.5GeV. 
Our results read: 

/_i/ + ,2# = -(0.23 ±0.02 ±0.01) GeV 3 , (134) 
f_if + ig' = -(0.41 ± 0.03 ± 0.04) GeV 5 (135) 

with the current (|5|), or 

/_ if+ ig' = (0.41 ± 0.03 ± 0.01) GeV 4 (136) 



with the current 



0.038 ±0.005 ±0.006 5 
f-,$f+.l9u= — GeV (137) 



* s u \ i 0.21 ±0.04 ±0.02 n 3 

/_ i/ + 3 [{r lt _ i + r 1>+ s)g + g ld ] = GeV 3 , (138) 

z z i i ilT'Q 

/-, 4 / + ,|[(r„,_,. +r 2 , + ,}) 9 "+ S f] = 0,17± °;° o 1±0 -° 2 GeV , (139) 
L,f + ^ + r^tf + S J] = ° 17±0 ° 1±0 02 GeV 3 . (140) 

11 i i m Q 

The central values conrespond to T = 1.3GeV and c<j c = 3.0GeV. The first and second 
error refers to the variation with T and uj c respectively. The inherent uncertainties due to 
the method of QCD sum rules and the pion wave functions are not included here. 

With the central values of f's in (|I~6|)-(|l8") and the values of r's in section |3], we arrive 

at: 

g = -(4.83 ±0.4 ±0.2) GeV" 2 , (141) 
g = -(4.43 ±0.33 ±0.44) (142) 



with the current (|6]), or 
with the current (01), 



g' = -(4.1 ±0.3 ±0.3) (143) 



0.77 ±0.19 ±0.19 
9u = , 144 

m Q 
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9u = - 2 -^^ GeV- 2 , (145) 

VTiQ 

= GeV _ 2 

rriQ 

= _u± M GeV _ 2 

We have not included the errors due to the uncertainty of f's in (|141|) - fll4"7|) , which is about 
20%. 

Therefore, with the pion wave function (p n (u) in ( |112| ) the D-wave coupling constants 
are g d (B h B*) = -(4.83 + J^) = -5.27GeV~ 2 , g d (B*,B) = -(3.94 + = -4.27GeV" 2 , 

^(S*, £*) = -(3.94 + ^i) = -4.34GeV- 2 with m b = 4.8GeV. Strong cancellation ocurrs 
for 0(1/i71q) corrections to the coupling constants from the operator /C. At the phenomeno- 
logical side the factor (r x _ i + r i,+,|) is somewhat large. Luckily the 0(l/mg) correction 
on the theoretical side has the same sign as (r 1 _ i + r 1+ -i)g. Moreover their magnitude 
is close to each other. Their contribution cancels largely. As the result the 0(1/771q) cor- 
rection is about of 9% for the B\ and 8% for B\ meson respectively. As mentioned above, 
since (|101|) , (|107|) and (|108|) depend on the second derivative of u^ip^iu) at Uq = |, the 
OiXjvriQ) correction to the D-wave coupling constants has a weak dependence on the the 
detailed shape of (p n (u) around uq = \. For example, if we use the pion wave function 
(p n {u) in ( |110| ), the C(l/mg) correction to the D-wave coupling constants is about 4% for 
the Bi and 3% for B\ meson respectively. The OiXjmQ) correction has the same sign as 
the leading order coupling constant g using all the five forms of (p w (u) in ( jl(J9| )-( |TT3"| ). 



9 Pionic widths of mesons in the (l + ,2 + ) doublet and 
discussions 

From the above results, it is straightforward to derive the decay widths for Bi and B\. 
Summing over charged and neutral pion channels the decay widths are: 

r(2?! - SV) = ±-{ 9 2 + 2gg ld )\g\ 5 + ^-gf s \q\ , (148) 
T(B* 2 -> Bn) = -Lfo2 + 2g a gl)\q\ 5 , (149) 

ZVTY 

T(B* 2 -> B*n) = -^-{gl + 2g h g\)\qf . (150) 

407T 

The last term in ( |148j ) is formally of 0(1/i71q). However, due to the large space phase 
factor for the S-wave decay it is numerically not small. 

The masses of Bi,B% are not known experimentally. The decay width is sensitive to 
the Bi, B\ masses. If we use the pion wave function ( |112j ) in []14[ |, the decay widths are: 

r < B - - bv > - ™ x (») 5 + 12 x (<4) 2 yy Mev • < i5i > 
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r (B ;^ B .) = 5.4x( i J|^) 5 MeV, (152) 

r ( B 2 «^ B V)=4.9x( i J^) 5 MeV, (153) 

where the decay momentum |g| is in unit of MeV respectively. 

If we assume the heavy quark expansion is also good for the charm system, the decay 
widths of Di and D\ are: 




r(Z>i -»• D*tt) = 6.2 + 9.7 x [ 7 ^rz\ MeV, (154) 



T(D* 2 -> Dti) = 13.7 MeV , (155) 

T(D* -> L>*tt) = 6.1 MeV , (156) 

where we have used the decay momenta |gj = 355, 503 and 387 Mev respectively for the 
above three processes listed in [p4| . Experimentally, the total decay widths of D\ and D2 
are 18.9l£;jjMeV and 23^MeV respectively. 

In summary, in this work we have improved previous results about the kinetic energy 
K, and chromomagnetic splitting £ of the excited heavy mesons in HQET. We have carried 
out a systematic expansion up to the order 0(1/ tuq) for the pionic decay amplitutes of 
Bi and B\ in the framework of HQET and light-cone sum rules. The mixing of the two 
1 + states has been properly taken into account. The 0{l/mq) correction to the decay 
amplitutes of B\ and B^ derived from the light-cone sum rule is small. When applied 
to D\ and D\ the central theoretical values are in good agreement with the experimental 
data, though the theoretical uncertainty is sizable for the D system. 
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Figure Captions 



Fig. 1. The sum rules for Gk as a functions of the Borel parameter T for 
different values of the continuum threshold u c . From bottom to top the curves 
correspond to uj c = 3.2, 3.0, 2.8 GeV respectively. The unit of T, Gk is GeV. 

Fig. 2. The sum rules for Gs with tu c = 3.2, 3.0, 2.8 GeV with the curves from 
top to bottom. 

Fig. 3. The sum rules for Gi a If , i with u c = 3.2,3.0,2.8 GeV. 

2>2 ~r, 2 

Fig. 4. The sum rules for /_ i / , zg with uj c = 3.2, 3.0, 2.8 GeV with the pion 

2 ti 2 

wave function in ( |112| ) in [14|. 

Fig. 5. The sum rules for /_ if, ig 1 with the current @ and lo c = 3.2, 3.0, 2.8 
GeV. 

Fig. 6. The sum rules for /' l f, id with the current (HI) uj c = 3.2,3.0,2.8 
GeV. 

Fig. 7. The sum rules for mqf_ i/ + zg\ s with u> c = 3.2,3.0,2.8 GeV. 



Fig. 8. The sum rules for m.Qf_ i/ + s y(r 1 _ i + r 1 + §_)g + gidj with oj c = 
3.2,3.0,2.8 GeV. 

Fig. 9. The sum rules for rn,g/_ i/ + 3 ^(r Q _ 1 + r 2 + 3)g a + gfj with w c = 
3.2, 3.0, 2.8 GeV. 

Fig. 10. The sum rules for rriQf_ if + 3 ur x _ 1 + r 2 + 3)9^ + g^j with 
w c = 3.2, 3.0, 2.8 GeV. 
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